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Abstract 

This article derives several properties of the Riesz distributions, such as their cor- 
responding Bartlett decompositions, the inverse Riesz distributions and the distri- 
bution of the generalised variance for real normed division algebras. In addition, 
introduce a kind of generalised beta distribution termed beta-Riesz distribution for 
real normed division algebras. Two versions of this distributions are proposed and 
some properties are studied. 

1 Introduction 

It is imminent the important role played by Wishart and beta distributions type I and II 
in the context of multivariate statistics. In particular, the relationship between these two 
distributio ns to obtain the beta dist ribution in terms of the distributi o n of t wo Wishart 
matrices. Faraut and Koranvi ( 1994 ). and subsequently lHassairi et al. ( 2005 ). propose a 



beta-Riesz distribution, which contains as a special case to the beta distribution obtained in 
terms of the distribution Wishart, which shall be named beta- Wishart, all this subjects in 
the context of simple Euclidean Jordan algebras. Such beta-Riesz distribution is obtained 
analog ously to the beta- W ishart distribution, but starting with a Riesz distribution. Re- 
cently, iDi'az-Garci'al (|2012t ) proposes two versions of the Riesz distribution for real normed 



division algebras. 

Based in these last two versions of the Riesz distributions, it is possible to obtain two 
versions of the beta-Riesz distributions, which by analogy with the beta- Wishart distribu- 
tions are termed beta-Riesz type I. As in classical beta- Wishart distribution, in addition it 
is feasible to propose two version for the beta-Riesz distribution type II. Each of the two 
versions for each beta-Riesz distributions of type I and II, (both versions for each) contain 
as particular cases to beta- Wishart distribution type I and beta- Wishart distribution type 
II, respectively. 

This article studies two versions for beta-Riesz distributions type I and II for real normed 
division algebras. Section [2] reviews some definitions and notation on real normed division 
algebras. And also, introduces other mathematical tools as two definitions of the generalised 
gamma function on symmetric cones, two Jacobians with respect to Lebesgue measure and 
some integral results for real normed division algebras. Section[3]proposes diverse properties 
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of two versions of the Riesz distributions as their Bartlett decompositions, inverse Riesz 
distributions and the distribution of the generalized variance. Section [4] introduces two 
generalised beta functions and, in terms of these, two beta-Riesz distributions of type I and 
II are obtained for real normed division algebras. Also, the relationship between the Riesz 
distributions and the beta-Riesz distributions are studied. This section concludes studying 
the eigenvalues distributions of beta-Riesz distributions type I and II in their two versions 
for real normed division algebras. 



2 Preliminary results 

A detailed discussion of real normed division algebras may be found in iBaed (|2002l) and 



Ebbinghaus et al. ( 1990h . For convenience, we shall introduce some notation, although in 



general we adhere to standard notation forms. 

For our purposes: Let F be a field. An algebra 21 over F is a pair (21; m), where 21 is a 
finite- dimensional vector space over F and multiplication m : 21 x 2t — > A is an F-bilinear 
map; that is, for all A € F, x, y, z € 21, 

m(x, Xy + z) — Xm(x\ y) + m(x; z) 
m(Xx + y; z) — Xm(x; z) + m(y; z). 

Two algebras (21; m) and (£;n) over F are said to be isomorphic if there is an invertible 
map cj) '■ 21 — > <£ such that for all x, y G 21, 

cj)(m(x, y)) = n(cj)(x), <p(y)). 

By simplicity, we write m(x; y) = xy for all x, y £ 21. 
Let 21 be an algebra over F. Then 21 is said to be 

1. alternative if x(xy) = (xx)y and x(yy) — [xy)y for all x, y £ 21, 

2. associative if x{yz) = (xy)z for all x,y,z £ 21, 

3. commutative if xy = yx for all x, y € 21, and 

4. unital if there is a 1 € 21 such that xl = x = lx for all x £ 21. 

If 21 is unital, then the identity 1 is uniquely determined. 

An algebra 21 over F is said to be a division algebra if 21 is nonzero and xy = Oa => ir = Oa 
or y = Oa for all x, y € 21. 

The term "division algebra" , comes from the following proposition, which shows that, in 
such an algebra, left and right division can be unambiguously performed. 

Let 21 be an algebra over F. Then 21 is a division algebra if, and only if, 21 is nonzero 
and for all a, b £ 21, with b ^ Oa, the equations bx = a and yb — a have unique solutions 
x,y G 21. 

In the sequel we assume F = 5i and consider classes of division algebras over or " 'real 
division algebras" for short. 

We introduce the algebras of real numbers 3?, complex numbers £, quaternions and 
octonions D. Then, if 21 is an alternative real division algebra, then 21 is isomorphic to 5i, 
£, .f) or D. 

Let 21 be a real division algebra with identity 1. Then 2t is said to be normed if there is 
an inner product (•, •) on 21 such that 

(xy, xy) = (x, x)(y, y) for all x, y £ 21. 
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If 21 is a real normed division algebra, then 21 is isomorphic 5ft, £, or O. 

There are exactly four normed division algeb r as: re al numbers (5ft), complex numbers 
(<£), quaternions ($)) and octonions (£)), see Baezj ( 2002 ). We take into account that should 
be taken into account, 5ft, £, 9) and D are the only normed division algebras; furthermore, 
they are the only alternative division algebras. 

Let 21 be a division algebra over the real numbers. Then 21 has dimension either 1, 2, 4 
or 8. In other bran c hes o f mat hemat i cs, th e parameters a = 2/(3 and t — /3/4 are used, see 



Edelman and Rao I ([20051 ) and iKabd (| 1984T) , respectively. 



Finally, observe that 
5ft is a real commutative associative normed division algebras, 
£ is a commutative associative normed division algebras, 
is an associative normed division algebras, 
D is an alternative normed division algebras. 
Let ££j n be the set of all to x n matrices of rank to < n over 21 with to distinct positive 
singular values, where 21 denotes a real finite- dimensional normed division algebra. Let 
2l mxn be the set of all to x n matrices over 21. The dimension of 21™ xn over 5ft is (3mn. Let 
A e 2l mx ™, then A* = A T denotes the usual conjugate transpose. 

Table Q] sets out the equivalence between the same concepts in the four normed division 
algebras. 



Table 1: Notation 



Real 



Complex 



Quaternion 



Octonion 



Generic 
notation 



Semi-orthogonal Semi-unitary Semi-symplcctic 



Orthogonal 
Symmetric 



Unitary 
Hcrmitian 



Symplcctic 

Quaternion 
hcrmitian 



Scmi-cxccptional 
type 
Exceptional 
type 
Octonion 
hcrmitian 



el 



,; ' m such that S = S*. In addition, 
Thus, consist of all matrices 



D = diag(di, . . .,d m ). 
71 such that ti4 = 



We denote by the real vector space of all S € 21 
let be the cone of positive definite matrices S € 21 
S = X*X, with X e £f nj „; then ^ m is an open subset of 6^. 

Let be the diagonal subgroup of C^ nm consisting of all D 6 2l m>< 

Let Ty(m) be the subgroup of all upper triangular matrices T 6 21'' 
for 1 < i < j < to. 

For any matrix X € 21™ xm , c?X denotes the matrix of differentials (dxij). F inally, we 
define the measure or volume element (dX) wh e n X € 2l mxn , 6^, or „, see Dumitriu 
(|2002h and D iaz-Garcia and Gutierrez-Jaimej (|201 ll ) . 

If X G 2l mxrl then (dX) (the Lebesgue measure in 2t mx,i ) denotes the exterior product 
of the j3nm functionally independent variables 



(dX) = A A dx^ where 



dxi 



f\ dx 



(fe) 



i=l j=l 



fc=l 



If S G (or S € T^r(m) with tjj > 0, i = 1, . . . , to) then (dS) (the Lebesgue measure in 
6^ or in T^(m)) denotes the exterior product of the exterior product of the m(m~ l)/3/2+m 
functionally independent variables, 



(dS) - /\ ds u /\ f\ ds 



(k) 



i<j k=l 
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Observe, that for the Lebesgue measure (dS) defined thus, it is required that S £ tymi that 
is, S must be a non singular Hermitian matrix (Hcrmitian definite positive matrix). 

If A e then (dA) (the Legesgue measure in S^J denotes the exterior product of the 
/3m functionally independent variables 

n 

(dA) = /\/\dAf ) . 

i=l fe=l 

If Hi e V£ >n then 

(HtdH,)= A A h J dh - 

i— 1 i+1 

where H = (H^H^)* = (hi, . . . , h m |h TO+ i, . . . , h„)* € &P(n). It can be proved that this 
differential form does not depend on the choice of the H2 matrix. When n = 1; 1 defines 
the unit sphere in 2l m . This is, of course, an (to — l)/3- dimensional surface in 2l m . When 
n — to and denoting Hi by H, (HcJH*) is termed the Haar measure on il^(m). 
The surface area or volume of the Stiefel manifold n is 

Vol«„) = / (HidH!) = (1) 

JH ie v£,„ [n/3/2] 

where T^Ja] denotes the multivariate Gamma function for the space 6^1 . This can be 
obtained as a particular case of the generalised gamma function of weight k for the space 
(5^ with k — (ki, &2> • • • i fcm), fci > > • • • > fc m > 0, taking At = (0, 0, . . . , 0) and which 
for Re(a) > (to — l)/3/2 — k m is defined by, see lGross and Richards! (Il987l) . 

I*[a,«] = [ ctr{-A}|Ar-( m - 1 )^ 2 - 1 g K (A)(dA) (2) 

m 

= 7r ™("»-i)/3/4j|r[ a + / Cj _( i _ 1)^/2] 

i=l 

= Mfr^[a], (3) 

where etr(-) = exp(tr(-)), | • | denotes the determinant, and for A € 

m—l 

q K (A) = |A m | fe ™ JT \A^- k ^ (4) 

i=l 

with An = (a r s), r, s = 1, 2, . . . ,p, p = 1, 2, . . . , to is termed the highest weight vector, see 
Gross and Richards! (|1987). Also, 



r£[a] = / etrl-AIIAI ^" 1 ^' 3 / 2 - 1 ^) 

m 

= n m(,m-i)0/iY[T[a-(i- 1)13/2], 



1=1 



and Re(a) > (to - l)/3/2. 

In other branches of mathematics the highest weight vector q K (A) is also terme d the 
generalised power o f A and is denoted as A K (A), see iFaraut and Koranvi (1994) and 



Hassairi and Laimil (|200ll ). 



Additional properties of q K (A), which are immediate consequences of the definition of 
q K (A) and the following property 1, are: 
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1. if A i, ... , A m , are the eigenvalues of A, then 

m 

9,(A)=n^- (5) 

i=i 

2. 

g K (A- 1 )= <z - 1 (A)= (Z _ K (A), (6) 

3. if k = (p, . . . ,p), then 

?«(A) = |A|", (7) 
in particular if p = 0, then q K (A) = 1. 

4. if r = (ti,*2> • • • ,*m), *i > h > ■ ■ ■ > t m > 0, then 

q K+T (A)=q K (A)q T (A), (8) 
in particular if r = (jp,p, . . . ,p), then 

fe+r (A) = q K+p (A) = \ A\ p q K (A). (9) 

5. Finally, for B e 2l mxm in such a manner that C = B*B € 6^, 

g K (BAB*) = <7 K (C)g K (A) (10) 

and 

q K (B- 1 AB*- 1 ) = (q K (C))- 1 q K (A). (11) 

Remark 2.1. Let P(&^ n ) denote the algebra of all polynomial functions on 6^, and 
Vki&m) the subspace of homogeneous polynomials of degree k and let V K (&^ 1 ) be an irre- 
ducible subspace of V{&^ n ) such that 

Note t hat q K is a homogeneous polynomial of degree k, moreover q K g V K {&^ n ), see Gross and Richards] 
(|l987h . 

In ([3|), [a]£ denotes the generalised Pochhammer symbol of weight k, defined as 

m 

[at = IKa-ft-l)^ 

i=l 

m 

^mCm-l)^ 7J r [ a + fc . _ (j _ l)p/2] 



rf„H 

r m[ Q ^] 

rt[a] ' 

where Re(a) > (m — l)/3/2 — fc m and 

(a)! = a(a + 1) • • • (a + i — 1), 
is the standard Pochhammer symbol. 
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An altern ative definition of the generalised gamma function of weight k is proposed by 

(12) 



Khatril ([19661 ). which is defined as 

= 7r »n("»-iW4jjr[a-jfe i -(m-i)j9/2] 

i=l 

(-i) fc rf>] 

[-o + (m - 1)^/2 + 



(13) 



where Re (a) > (to - 1)0/2 + fc j. 

Similarly, from iHera ( 1955 , p. 480) the multivariate beta function for the space S^, can 
be defined as 

B&M = / |S| Q -( m - 1 ^/ 2 - 1 |I m -S| b -( m - 1 )' 3 / 2 - 1 (dS) 
io<s<i m 



ri[a + b] ' 



(14) 



where R = (I— S) 1 — I, Re(a) > (to— 1)0/2 and Re(6) > (m— 1)0/2, see Dfaz-Garcia and Gutierrez- Jaimezl 
(|2009bh . 

Now, we show three Jacobians in terms of t he parameter, which are proposed as ex - 
tensions of real, complex or quaternion cases, see lDi'az-Garci'a and Gutierrez- Jaimezl (|201ll ). 

Lemma 2.1. Let X and Y 6 matrices of functionally independent variables, and let 
Y = AX A* + C, where A € ££, „, and C 6 Q3£, are matrices of constants. Then 

' i It, , i II 'ill •> 



(dY) = |A*A|^( m - 1 )/ 2+1 (dX). 



(15) 



Lemma 2.2 (Cholesky's decomposition). Let S € and T g Ty-(m) with tu > 0, 
z = 1,2,..., to. TTien 



(dS) = 



2 m Y[tii m ~ i)+1 ( dT ) if s = T * T ; 

i=l 
m 

2 m JJtf l (l ~ 1)+1 (dT) i/S = TT*. 



(16) 



Lemma 2.3. Let S € tym- Then, ignoring the sign, if Y = S 1 

(dY) = |ST' 3 ( m - 1 )- 2 (dS). 



(17) 



Next is stated a general result, that is useful in a variety of situations, which enable 
us to transfor m the density function of a matrix X € ffi ff, to the density function of its 
eigenvalues, see Dfaz-Garcfa and Gutierrez- Jaimez (l2009bl) . 

Lemma 2.4. Let X 6 be a random matrix with a density function / X (X). XTien the 
joint density function of the eigenvalues Ai, . . . , A m o/ X is 



r m 2 /3/2+f> 



r& [to/3/2] 



m „ 

J](A. 1 -A J )' 3 / /(HLH*)(dH) 



(18) 
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where L = diag(Ai, . . . , A m ), Ai > • • • > A m > 0, (dH) is the normalised Haar 



measure ana 



Q = 



0, = 1; 

-to, (3 = 2; 

-2m, = 4; 

-4m, = 8. 



Finally, let's recall the multidimensional convolution theorem in terms of the Laplace 
transform. For this purpose, let's use the complexification = &^ + i&t of &t- That 
is, consist of all matrices X € (g<)™x™ of the form Z = X + iY, with X, Y G 

It comes to X = Re(Z) and Y = Im(Z) as the rea/ and imaginary parts of Z, respectively. 
The generalised ri ght half-plane = y£, + z<S£, in 6^j C consists of all Z € 6jJ; £ such that 



Re(Z) G see ([Gross and Richards! . Il987. p. 801). 



Definition 2.1. If /(X) is a function of X € ^3^, the Laplace transform of /(X) is defined 
to be 

g(T) = ( etr{-XZ}/(X)(dX). (19) 

where Z G 

Lemma 2.5. // <?i(Z) and 32(2) are f/ie respective Laplace transforms of the densities 
/ X (X) and ff Y (Y) £/ien £/ie product g\{Ti)gi{Ti) is the Laplace transform of the convolution 

/x( X ) *5yOO; where 

^ E (S) = / x (X)*p Y (Y)= / / x (T) ffv (S-T)(dT), (20) 

•/0<T<H 

urafft S = X + Y and T = X. 

3 Riesz distributions 



This section shows two versions of the Riesz distributions (jDiaz-Garcfal |2012) and the study 



of their Bartlett decompositions. Also, inverse Riesz distributions are obtained. 
Definition 3.1. Let S G 3?^ and k = {k\,ki, . . . , k m ), ki > k% > • • ■ > k m > 0. 

1. Then it is said that X has a Riesz distribution of type I if its density function is 

for Re(a) > (m — 1)0/2 — k m ; denoting this fact as X ~ ^^(a, k, S). 

2. Then it is said that X has a Riesz distribution of type II if its density function is 

_ etr{-S- 1 X}|Xr-(" l - 1 )' 3 / 2 - 1 g K (X- 1 )(dX) (22) 

rf„[a,- K ]|S| a 

for Re(a) > (m - 1)0/2 + fcr; denoting this fact as X — Ui^i II (a, n, E). 

Theorem 3.1. Let £ G 3?^ and n = (fci, &2, . . . , fc m ), k\ > &2 > • • • > k m > 0. j4nd fe< 
Y = X- X . 
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1. Then i/X has a Riesz distribution of type I the density of Y is 

aom+^^j ki 

/ r— — etr{-E- 1 Y- 1 }lYl-( a +C">- 1 )^ a + 1 ) gw (Y- 1 )(dY) (23) 
ri[a,K]\-S\ a q K CE) 

for Re(a) > (m — l)/3/2 — fe m and is termed as inverse Riesz distribution of type I. 

2. Then if X has a Riesz distribution of type II the density of Y is 

etri-S-^-HlYl-^+^-^^+^^fY)^) (24) 

r^[a,- K ]|S| a 

for Re(a) > (m — 1)13/2 + k\ and it said that Y has a inverse Riesz distribution of 
type II. 

Proof. It is immediately noted that (dX) = | Y| -■ s C T ™- 1 >- 2 (dTY) and from $ZT$ and (J22). □ 

□ 

Observe that, if k = (0, 0, . . . , 0) in two densities in Definition 13.11 and Theorem 13. II the 
matrix multivariate gamma and inverse gamma distributions are obtained. As consequence, 
in this last case if £ = 2S and a = /3n/2, the Wishart and inverse Wishart distributions 
are obtained, too. 

Theorem 3.2. Let T e 1^(m) with tu > 0, % = 1, 2, . . . , m and define X = T*T. 

1. //X has a Riesz distribution of type I, \21\) , with X = I m , then the elements Uj (1 < 
i < j < m) of T are all independent. Furthermore, t^ ~ G (a + ki — (i — l)/3/2, 1) 
and V2tij ~ A/f (0, 1) (1 < i < j < m). 

2. //X has a Riesz distribution of type II, \22\i , with S = l m , then the elements Uj (1 < 
i < j < m) of T are all independent. Moreover, ~ (a — ki — (i — l)/3/2, 1) and 
V2tij ~ A/f (0, 1) (1 < i < j < m). 

Where x ~ G^{a, ci) denotes a gamma distribution with parameters a and a and y ~ A/"-f (0, 1) 
denotes a random variable with standard normal distribution for real normed division alge- 
bras. Moreover, their respective densities are 

Q\x : a,a) = 1 expl-Z^/aK- 1 ^), 
(a//3) Q I>] 

and 

<(2/ : 0, 1) = (2 , ^ exp{-/3y 2 /2}(^) 

where Re(a) > and a G 3>i , see lDzaz-Gama and Gutierrez- Jdimez It 201 A ). 

Proof. This is given for the case of Riesz distribution type I. The proof for Riesz distribution 
type II is the same thing. The density of X is 

P , - etr{-X}|X| a -("'- 1 )' 3 / 2 - 1 g K (X)(dX). (25) 
Tfn[a, k] 



s 



Since X = T*T we have 

m 

trX = trT*T = ^i|, 

i<j 

m 

|X| = |T*T| = |T| 2 =n4, 

in — 1 m 

q K (X) = q K (T*T) = |T*T| fc ™ [] |T|Tj| fei ~ fei+1 = JJ^**, 

i=i i=i 

and by Theorem 12.21 noting that dt\ = 2tadtu, then 

m / 

(dx) = 2 ™n4 (w ^ )+i A** 



= n(^) M2 (A^)A(A^ 

i=l \i=l / yi<j 

Substituting this expression in (p?5j) and using ([5]) we find that the joint density of the 
tij (1 < i < j < m) can be written as 

m aa+ki-{i-l)P/2 .. . 

n tH-o-M ^ (lJ) Wi) 

m ^ 

x II ( 7r / j8 )/3/2 expi-^-X^j), 

only observe that 



g a +ki-(i-l)P/2 m i 



□ 



□ 



In analogy to generalised variance for Wishart case, the following result gives the distri- 
bution of |X| when X has a Riesz distribution type I or type II. 

Theorem 3.3. Let v = |X|/|S|. Then 

1. z/X has a Riesz distribution of type I, h21}) , the density of v is 



(tl:a + h-(i-l)p/2,l). 

i=i 

2. if X has a Riesz distribution of type II, h22]) , the density of v is 

m 

[]^ {t 2 u :a-h-(i-l)p/2,l). 
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where t 2 it i = l,...,m, are independent random variables. 
Proof. This is immediately from Theorem 13.21 noting that if 

B = E" 1/2 XS" 1/2 = T*T, 
with T 6 %^(m) and ijj > 0, i = 1, 2, . . . , m, then 



|B| = jj^ = |X|/|E|=t,. □ 



□ 

4 Generalised beta distributions: Beta-Riesz distribu- 
tions. 

This section defines several versions for the beta functions and their relation with the gamma 
functions type I and II. In these terms, the beta-Riesz distributions type I and II are defined. 
Finally, diverse properties are studied. 

4.1 Generalised c-beta function 

A generalised of multivariate beta function for the cone ^3^, denoted as B m [a, k; b, r], can 
be defined as 

f \S\ a ^ m - 1 ^/ 2 - 1 qK (S)\I m -S\ b -^- 1 ^/ 2 - 1 q T (I m -S)(dS) (26) 

Jo<s<i m 

where n = (fci, k 2 , . . . , k m ), fci > fe > ■ • • > k m > 0, r = {h,t 2 , . ■ .,t m ), h > t 2 > ■ ■ ■ > 
t m > 0, Re(a) > (m - 1 )0/2 - k m and Re (6) > (m - l)/3/2 - t m . This is defined by 



Faraut and Koranvil ( 19941 . p. 130) for Euclidean simple Jordan algebras. In the context of 



multivariate analysis, this generalised beta function can be termed generalised c-beta function 
type I, as analogy to the correspondence case of matrix multivariate beta distribution, and 
using the term c-beta as abbreviation of classical-beta. In the next theorem we introduce 
the generalised c-beta function type II and its relation with the generalised gamma function. 

Theorem 4.1. The generalised c-beta function type I can be expressed as 

| R |a- (m -l)/V2-l gK(R) | Im + R \-(a+b) q -l r{Im + R)(dR) 



L 



_ r m [a,K]T^[b,r} 



T m \a + b, k + t] 



where K= {k\, k 2 , ■ ■ ■ , k m ), k\ > fc 2 > • • • > k m > 0, r = (ti,t 2 , ■ . ■ ,t m ), h > t 2 > ■ ■ ■ > 
t m > 0, Re(a) > (m — 1)0/2 — k m and Re{b) > (m — 1)0/2 — t m . The integral expression is 
termed generalised c-beta function type II. 

Proof. Let R = (I S)" 1 - I then by Lemma H (dS) = (I + R)~' 3(m - 1) - 2 (d(R)). The 
desired result is obtained making the change of variable on (|26| noting that S = (I m + 
R)- 1 / 2 R(I m + R)" 1/2 , then 

|S| = IRH^ + RI" 1 , 
IL-S^IL+RI- 1 , 
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and by ^ and (ITTt we have 

q K ((I m + R)~ 1//2 R(I TO + R) _1//2 )g T ((Im + R)" 1 ) = q K (R)q-l T (I m + R). 

For the expression in terms of generalised gamma function, let B = H 1 / 2 SH 1 / 2 in (|26|) . such 
that (H 1 / 2 ) 2 = S. Then (dS) = |H|-( m " 1 )' 3 / 2 - 1 (dB), and 

^[a, K ;6,r]|Hr+ b -(™- 1 W 2 - 1 (7K+T (H) 

= f |B| Q -( m - 1 ) /3 / 2 - 1 g K (B)|H-B| b -( m - 1 ^/ 2 - 1 g T (H-B)(rfB). 

JO 

Taking Laplace transform of both size, by ([2~Tj) . the left size is 

/ B£[a, k; b, r] etr{-HZ}|Hr+ fc -(" l - 1 )^ 2 - 1 (ZK+T (H)(rfH) 

= ^[ a , k; b, r]lt [a + 6; k + t]\Z\-^ q K+r {Z), 

and applying Lemma 1231 <7i(Z) is 

/ etr{-SZ}|H|«-(— l W 2 -\ K (3)(dS) = I* [a;K]|Z|-" ?(t (Z), 
and <?2(Z) is given by 

/ etr{-HZ}|H - B^-D^-l^g _ B)(dS) = r £ [6; r ]| Z |-V(Z). 



Thus, equally 



1 Ki a + o, k + t 



□ 



4.2 Generalised A>beta function 



Alternatively, a generalised of multivariate beta function for the cone denoted as 

B^[a, — k; b, — r], can be defined as 



/ | Sr -(— D/3/2-l (?K(S -l ) | Im _ S | fe -( OT -l)^/2-l ?r((Im _ S) -l )(dS) 

^0<S<I m 



(27) 



where k = (fei, fc 2 , . . . , k m ), ki > k 2 > • • • > k m > 0, r = (ii, t 2 , ■ ■ -,t m ), h >t 2 > ■■■ > 
t m > 0, Rc(a) > (to - l)/3/2 + fci and Rc(6) > (to - l)/3/2 + Again, in the context 
of multivariate analysis, this generalised fc-beta function can be termed generalised k-beta 
function type I, as an analogy to the corresponding case of matrix multivariate beta distri- 
bution and using the term k-beta as abbreviation of Khatri-beta. Next theorem introduces 
the generalis ed k-be t a func tion type II and its relation with the generalised gamma function 
proposed by iKhatril ([1966). 

Theorem 4.2. The generalised k-beta function type II can be expressed as 

f |R|a-(m- lW2-l g -l( R) | Im + K\-^q K+T (l rn + R)(dR) 
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_ Ti[a,- K }ri[b,-r} ^ 
T m [a + b, —k — t] 

where K — (fci, k 2 , ■ ■ ■ , k m ), k\ > k% > ■ ■ ■ > k m > 0, r = (ti,*2, • • • ,t m )> h > h > • • • > 
t m > 0, Re(a) > (m — l)/3/2 + k\ and Re(b) > (m — l)/3/2 + t\. The integral expression is 
termed generalised fc-beta function type II. 

Proof. The proof is analogous to the given for Theorem 14.11 □ □ 

Observe that if k = (0, . . . , 0) and r = (0, . . . , 0) in (J26J) , Theorem |4~T| ([27]) and Theorem 
1421 the beta function (fT4l) is obtained. 



4.3 c-beta-Riesz and fc-beta-Riesz distributions 

As an immediate consequence of the results of the previous section, next the c-beta-Riesz 
and fc-beta-Riesz distributions types I and II are defined. 

Definition 4.1. Let k = (k%, . . . , k m ), k\ > k 2 > ■ ■ ■ > k m > and r = (t\,t2, ■ ■ ■ , t m ), 
h>h>---> t m > 0. 

1. Then it said that S has a c-beta-Riesz distribution of type I if its density function is 

isia-Cm-l^-^g^ _ g 1 6- ( m - 1)/3/2- 1 g ^ (I ^ _ (2g) 

B m [a, k;o,t] 

where < S < I m and Re(a) > (to - l)(3/2 - k m and Re(6) > (m — l)/3/2 - t m . 

2. Then it said that R has a c-beta-Riesz distribution of type II if its density function is 

1 iRl'-^Wa-ig^)^ + R|-(«+f) g -i T (I m + R)(dR), (29) 
Dm [a, k; 6, rj 

where R € *pf„ and Re(a) > (m - l)/3/2 - fc m and Re(6) > (m - l)/3/2 - i m . 
Similarly we have 

Definition 4.2. Let k = (k%, k 2 , . . . , fc m ), fci > k 2 > ■ ■ ■ > k m > and r = (£i,t2j ■ • ■ j im), 
ii > t 2 > ■ ■ ■ > t m > 0. 

1. Then it said that S has a k-beta-Riesz distribution of type I if its density function is 

1 JslMm-U/J^-lg-l^i^ _ S |6-( m -l)^/2-l g -l (Im _ s)(ds)j (3Q) 

Bm[o, -k; 6, -r] 

where < S < I m and Re(a) > (m - l)/3/2 + fci and Re(6) > (m - l)/3/2 + t x . 

2. Then it said that R has a k-beta-Riesz distribution of type II if its density function is 

„ Pr 1 u -iRr-^-^-^CRJIIn. + Rl-^g^+r^+R)^), (31) 

Dm [a, -re; 6, -rj 

where Re^ and Re(a) > (m - l)/3/2 + fci and Re(6) > (to - l)/3/2 + h. 



Observe that the relationship between the densities (|28|) and (|29|). and between the 
densities (1501) and (pTTj) are easily obtained from the theorems 14.11 and 14.21 respectively. 
The following result state the relation between the Riesz and beta-Riesz distributions. 
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Theorem 4.3. Let Xi and X2 be independently distributed as Riesz distribution type I, 
such that Xi — $R^ J (a, K, I m ) and X 2 — 9i^j 7 (6, r, I m ) 7 i?e(a) > (to - l)/3/2 + fci and 
i?e(6) > (to - 1)/3/2 + *i. Let 

S = (Xi + x 2 )- 1 /2 X2 (x 1 + X2)- 1 / 2 ', 

w/iere (Xi + X 2 ) x / 2 is any nonsingular factorisation of (Xi + X 2 ), (Xi + X 2 ) = (Xi + 
X 2 ) 1 / 2 (Xi + X 2 ) 1 / 2 . Then S has a c-beta-Riesz distribution type I. 

Proof. The joint density of Xi and X 2 ) is given by 

^a+b)m+J2T=i( k i+ t i) 



r^[a,/t]r^[6,r] 



etr{-(Xi + X 2 )}|X 1 | Q -( m - 1 )^ 2 - 1 gK (X 1 ) 



x |X 2 | b -(™- 1 )' 3 / 2 - 1 g T (X 2 )(rfX 1 ) A (dX 2 ). 

Let Y = Xi + X 2 and Z = X 2 , then, (dXi) A (dX 2 ) = (dY) A (dZ). Then the joint density 
of Y and Z is given by 

§ ^ etr{-Y}|Y - Z| a -(™- 1 )' 3 / 2 - 1 o K (Y - Z) 

[a, K]r^,[6, r] 

x | Z | b -( m - 1 W 2 - 1 g r (Z)(cfY) A (dZ). 
Let S = Y- 1 / 2 ZY- 1 / 2 ' and W = Y, with Y = Y^'y 1 / 2 , then, 

(dY) A (dZ) = |W|"< ,n - 1 >/ 2+1 (dW) A (dS) 

Hence the joint density of S and W is 

i-^g ? etr{-W}|W| a+h -^ m - 1 )/ 2 - 1 o K+T (W)|I - S | a -( m - 1 )' 3 / 2 - 1 

[a,K]T p m [b,r] 

xq K (I - S)|S| b -( m - 1 )^ 2 - 1 o T (S)(dW) A (dS). 
Therefore, integrating with respect to W the desired result is obtained. □ □ 

Theorem 4.4. Let Xi and X 2 6e independently distributed as Riesz distribution type I, 
such that Xi ~ ^ / (a, K,I m ) and X 2 - fK^ 7 (6, t, I m ), i?e(a) > (to - l)/3/2 + fci and 
Re(b) > (to - 1)/3/2 + *i. Let 

r = x~ 1/2 x 2 x~ 1/2 ', 

where X/ is any nonsingular factorisation 0/X1, to i/ie sense 
S /las a c-beta-Riesz distribution type II. 

Proof. The joint density of Xi and X 2 ) is given by 

a(a+b)m+J2T=i( k i+ t i) 

UMliftlM ^-.x. + x,)) l x 1 r-.--)^-„(x 1 ) 

x|X 2 | b -( m - 1 )^ 2 - 1 g T (X 2 )(dX 1 ) A (dX 2 ). 
Let R = X~ 1/2 X 2 X~ 1/2 ' and W = X l7 with X x = Xj /2 'xJ /2 , then, 

(dXi) A (dX 2 ) = (Wl^-^^+^dW) A (dR) 
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Hence the joint density of S and W is 

P ~a tL 1 etr{-W(I m + RjJlw^MWu^fw) 

Lm[a, K\lm[0, T\ 

x | R |6-(m-l)/3/2-l (?r(R)(dw) A (dR ^ 

Then, integrating with respect to W the density of R is obtained. □ □ 



The following theorems 14.51 and 14.61 contain versions for /c-beta-Riesz distributions of 
theorems 14.31 and 14.41 whose proofs are similar. 



Theorem 4.5. Let Xi and X 2 be independently distributed as Riesz distribution type II, 
such that Xi ~ £R^ /J (a, k, I m ) and X2 ~ £H(^ 77 (fe, r, I m ), Re(a) > (m — l)/3/2 + fci and 
i?e(fe) > (to - l)/3/2 + Let 

S = (Xi + X 2 )- 1 / 2 X 2 (Xi + Xa)- 1 / 2 ', 

where (Xi + X2) 1 / 2 is any nonsingular factorisation of (Xi + X2), (Xi + X 2 ) = (Xi + 
X 2 ) 1 / 2 (Xi + X 2 ) x ^ 2 . TTien S /ias a k-beta-Riesz distribution type I. 

Theorem 4.6. Lei Xi and X 2 fee independently distributed as Riesz distribution type I, 
such that Xi - ^ u {a,K,l m ) and X 2 - <K^ 77 (fe, t, I m ), Re{a) > (to - 1)^/2 + ki and 
Re(b) > (to - l)/3/2 + ti. Let 

— A.j A2A; , 

where X^ 2 is any nonsingular factorisation 0/X1, in t/ie sense Xi = X*^ 2 Xj^ 2 . T/ien S 
/ias a k-beta-Riesz distribution type II. 

4.4 Some properties of the c-beta-Riesz and &;-beta-Riesz distribu- 
tions 

This section derives the distributions of eigenvalues for c-beta-Riesz and fc-beta-Riesz dis- 
tributions type I and II. First remember that: 

Remark 4.1. If chj(M) denotes the i-th eigenvalue of the matrix M £ Cp^, then note that 

ifx.Yeq^, 

ch l (X- 1/2 YX- 1/2 ) = chi(X _1 Y) = ch^YX" 1 ). 
Furthermore, if A e 3?^ is a nonsingular matrix then 

ch,((A'XA)- 1 A'YA) = chi(X -1 Y) 

Then, it is found the joint distributions of the eigenvalues for random matrices c- and 
/c-beta-Riesz type I and II assuming that X 7^ I m in theorems I4.3M.6I 

Theorem 4.7. Let X £ <&^, n = (k±, fc 2 , . . . , k m ), ki > k 2 > ■ ■ ■ > k m > and t = 

(ii,*2,...,im), h>t 2 >--->t m >0. 

1. Let Ai, . . . , X m , Ai > • • • > A m > fee the eigenvalues of S. Then if S has a c-beta- 
Riesz distribution of type I, the joint density of Ai, . . . , A m is 

m 2 13 /2+o m m 

— — n(Ai - n K +k ^ [m ' i)m ^ 



T^[mP/2]B^[a,K;b,T] 



i<j i—l 
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i=l \i=l / 

where < Aj < 1, i = 1, . . . , m and Re(a) > (m — l)/3/2 — k m and Re(b) > {in — 
1)0/2 -t m . 

2. Let 8i, . . . ,5 m , 8i > ■ ■ ■ > S m > be the eigenvalues o/R. Then i/R has a c-beta-Riesz 
distribution of type II, the joint density of their eigenvalues is 

m 2 B/2+n m m 



m / rn \ 

n(i-^r (o+i,+fei+i<) ' 

i=l \i=l / 

where 8i > 0, i = 1, . . . , m and Re{a) > (rn — l)/3/2 — fc m and Re(b) > (m— 1 ) /3/ 2 — t „ 
g is defined in Lemma\ 



Proof. This is due to applying the Lemma 12.41 in (|2"51) and (|2T)1) , and taking into account the 
Remark |4~T1 and equations (0 and □ □ 

This section conclude establishing the Theorem 14.71 for the case of the fc-beta-Riesz 
distributions. 

Theorem 4.8. Let X <E 3?^, k = (ki, fe, • • • , fc m ), fci > > • • • > fc m > and t = 

{t\,t2, ■ ■ ■ , t m ), ti > t 2 > ■ ■ ■ > t m > 0. 

1. Let Ai, . . . , A m , Ai > • • • > A m > be the eigenvalues of S. Then if S has a k-beta- 
Riesz distribution of type I, the joint density of X\ , . . . , X m is 



TT(Ai - a^TTa 

r^[m/3/2]^[ a ,- K ;6,-r]^ V 



a-fc;-(m-l)/3/2-l 



jj (1 _ Ai) ^-(™-iW2-i /y dA \ 

i=l \i=l / 

where < A^ < 1, i = 1, . . . , m and Re(a) > (rn — l)/3/2 + fci and i?e(o) > (rn — 
1)0/2 + t t . 

I?. Lei . . . , 5 m , 5\ > ■ ■ ■ > S m > be the eigenvalues of R. Then if R /ras a k-beta- 
Riesz distribution of type II, the joint density of their eigenvalues is 



-T 



T p m [mj3/2]B p m [a,~K;b, 



jrr - 3 II(* " II sr^-^ 2 - 1 



i<j i=l 



i=l \i=l / 

where 6i > 0, % = 1, . . . , m and Re(a) > (rn — l)/3/2 + fci and Re{b) > (rn — l)/3/2 + ti. 

g is defined in Lemma \2.4\ 

Finally observe that if in all result of this section are taking k — (0, 0, . . . , 0) and 
t = (0,0,..., 0) the obtained results are the corresponding to matrix multivariate beta 
distributions of type I and II. 
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Conclusions 



Finally, note that the real dimension of real normed division algebras ca n be express ed as 
powers of 2, j3 = 2™ for n = 0, 1, 2, 3. On the other hand, as observed from iKabd (11984 . the 
results obtained in this work can be extended to hypercomplex cases; that is, for complex, 
bicomplex, biquaternion and bioctonion (or sedenionic) algebras, which of course are not 
division algebras (except the complex algebra). Also note, that hypercomplex algebras are 
obtained by replacing the real numbers with complex numbers in the construction of real 
normed division algebras. Thus, the results for hypercomplex al gebras are o btained by 
simply replacing /3 with 2/3 in our results. Alternatively, following iKabel (Il984h . it can be 
concluded that, results are true for '2™-ions', n = 0,1,2,3,4,5, emphasising that only for 
n = 0, 1, 2, 3 are the result algebras, in fact, real normed division algebras. 
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